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We present a derivation of two-point correlations of general tracers in the peak-background split 
(PBS) framework by way of a rigorous definition of the PBS argument. Our expressions only 
depend on connected matter correlators and "renormalizcd" bias parameters with clear physical 
interpretation, and are independent of any coarse-graining scale. This result should be contrasted 
with the naive expression derived from a local bias expansion of the tracer number density with 
respect to the matter density perturbation 5l coarse-grained on a scale Rl- In the latter case, the 
predicted tracer correlation function receives contributions of order (52) at each pcrturbative order 
n, whereas, in our formalism, these are absorbed in the PBS bias parameters at all orders. Further, 
this approach naturally predicts both a scale-dependent bias cx such as found for peaks of the 
density field, and the scale-dependent bias induced by primordial non-Gaussianity in the initial 
conditions. The only assumption made about the tracers is that their abundance at a given position 
depends solely on the matter distribution within a finite region around that position. 
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I. INTRODUCTION 

The clustering of tracers of the large-scale structure 
(LSS) in the Universe, such as galaxies, clusters, or the 
Lyman-Q! forest, is one of the most important probes of 
the origin and evolution of cosmological perturbations. 
On sufficiently large scales, correlations are weak and one 
should hope that a perturbative approach will allow us to 
cleanly connect observations with the predictions of cos- 
mological models (such as the standard cold dark matter 
scenario with adiabatic Gaussian initial conditions). In 
case of the matter density perturbations, these pertur- 
bative approaches arc well developed (sec [1] for a re- 
view). On the other hand, the formation of the trac- 
ers we actually observe necessarily involves highly non- 
linear, small-scale mechanisms which cannot be described 
perturbatively. In order to relate observations to theory, 
we thus need an effective description which involves un- 
known "bias parameters". These parameters in general 
need to be determined observationally, and we would like 
as many parameters as necessary to accurately describe 
tracer correlations down to some minimum scale, but not 
more, in order to retain the maximum amount of cosmo- 
logical information. Which bias parameters need to be 
included, and to what order, is still an open problem. 

The simplest and most well-known bias expansion is 
the local expansion in terms of density [2-4], 

5h,L{^) = Co + ci(5i(x) + |(5i(x) + . . . , (1) 

where c„ are the bias parameters, 5h,L is the fractional 
tracer density perturbation, and 5^ is the correspond- 
ing matter density perturbation. Here, both the tracer 



and the matter density field are understood as coarse- 
grained on some scale Rl, so that this expansion can 
be interpreted as a counts-in-cells relation, and the c„ 
as "scatter-plot bias parameters" . This perturbative de- 
scription is commonly assumed to be valid on large scales 
[5, 6], i.e. only if a\ — ((5|) ^ 1, although the series 
could actually converge even for ctl ^ 1 if the condition 
\cn+i/cn\ ~ const, is satisfied in the limit n — >■ 00; how- 
ever, the larger ctl, the more terms need to be included 
to obtain a converged expression. As can easily be seen, 
the correlation function predicted by Eq. (1) depends on 
the coarse-graining scale through the variance cr£ = {5\), 
and all higher moments of the density field. 

In contrast to counts-in-cells studies [3], where a spe- 
cific scale Rl is singled out, wc expect that no additional 
smoothing scale Rl should enter the calculation of corre- 
lation functions on a scale r (unless we directly link Rl 
to r, as done in [7]). Thus, we need to absorb the Rl- 
dependent pieces into renormalized bias parameters fojv, 
as proposed for the first time by McDonald [8] and tested 
against simulations in Jcong and Komatsu [9]. More re- 
cently, this approach has been pursued to higher order in 
the multipoint propagator framework [10] by Matsubara 
[11]. In accordance with renormalization theory, the ex- 
pression for tracer correlators in terms of the parameters 
6jv should be i?i,-independent, and convergent in terms 
of matter correlators at the separation at which we mea- 
sure tracer correlations (rather than zero-lag correlators 
such as (Sl)). In the process of renormalization, we have 
to introduce a new bias parameter 6jv at each order in 5l 
which needs to be determined from observations. Simi- 
lar conclusions hold when adding other, non-local quan- 
tities to the expansion Eq. (1), such as derivatives of d 
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or the tidal tensor. With in principle arbitrarily many 
free parameters, the clustering of LSS tracers loses all 
cosmological constraining power. Moreover, treating the 
bN as mere nuisance parameters (as in [8, 12]) precludes 
us from using any information contained in the 6jv on the 
parent halos and formation history of the tracers. 

There is thus strong motivation to try to associate 
physical meaning with the renormalized bias parameters. 
This would allow models of, say, galaxy formation to pro- 
vide sensible priors on the allowed range of the param- 
eters, or to predict connections between different bias 
parameters. For example, there is a well- motivated con- 
nection between the linear bias with respect to density 
of a tracer, and the parameter that quantifies the scale- 
dependent bias induced by primordial non-Gaussianity 
[13, 14]. In this case, the connection is crucial, as with- 
out it one would not be able to constrain primordial non- 
Gaussianity from the scale-dependent bias. 

Our goal in this paper is to explicitly derive the physi- 
cal meaning of the renormalized bias parameters b^, and 
their exact relation to observables such as the correlation 
function to any order. Physically, the renormalized bias 
parameters quantify the response of the mean abundance 
of tracers to a change in the background matter density p 
of the Universe (at fixed proper time since the Big Bang), 
i.e. 



This definition applies to tracers of any nature; in the 
considerably more restrictive case of universal mass func- 
tions, where the abundance of tracers depends only on 
'^c = Sc/cr{M), a fractional change D in the back- 
ground density is equivalent to a change in the parameter 
Sc ~> 5c — leading to the well-known peak-background 
split (PBS) bias parameters [15], Eq. (49) in Sec. II C. 
Even though Eq. (2) is much more general, we will still 
refer to the parameters b^ as "PBS biases" for conve- 
nience, although it is important to keep the distinction 
in mind. 

We will show that the bias parameters defined through 
Eq. (2), and its generalization to non-local biases, are the 
coefficients multiplying powers of the matter correlation 
function ^^(f') in the expansion of the tracer correlation 
function. Specficially, for a Gaussian density field we 
obtain 

We also show that the 6jv agree with the bias parame- 
ters identified through a direct calculation of the correla- 
tions of thresholded regions and peaks of the density field. 
However, it is important to note the difference in philos- 
ophy between each derivation: thresholded regions and 
peaks constitute "microscopic" models of tracers, where 
the relation between matter and tracer density is explic- 
itly specified on all scales (this is also the approach taken 



in Matsubara [11], although the description is kept gen- 
eral there). On the other hand, here we treat local bi- 
asing (and its generalization to non-local quantities) as 
effective description on sufficiently large scales, indepen- 
dently of the microscopic physics. In the language of field 
theory, the former approaches constitute specific "UV- 
complete" theories, while the approach presented here is 
an "effective theory" of biasing. Specifically, Rl serves 
the role of a UV cutoff whose precise value should not 
impact correlations on scales of observable interest. 

The exact relation between the parameters 6jv and 
tracer correlations provides a rigorous framework in 
which further modeling assumptions for any given tracer, 
for example from the excursion set, peak model, or halo 
occupation distribution, can be embedded both in or- 
der to tighten constraints on cosmological parameters, 
and in order to infer the physics of the formation of the 
tracers (as also pointed out in [16], the deviations from 
the peak-background split predictions found in [17] are 
due to the fact that the authors assumed a universal mass 
function of halos, not due to the inaccuracy of the peak- 
background split argument itself). Along the way, we also 
show that renormalization removes the zero-lag matter 
correlators from the expression for the tracer correlation 
function at all orders, as required. We also show that 
the same bias parameters describe both the tracer auto- 
and the cross-correlation with matter. The expression of 
tracer correlations in terms of renormalized bias param- 
eters and "no-zero-lag" matter correlators such as the 
correlation function is manifestly convergent as long as 
these matter correlators are small. This is in close anal- 
ogy to the expression for tracer correlations in terms of 
resummed bias parameters in [11]. Note that the treat- 
ment in [11] is in Fourier space, while we work in real 
space here, for which we believe that the physical as- 
sumptions and arguments are more clear. 

The renormalization approach proposed in [8] provides 
us with another extremely valuable tool: in describing 
the tracer density in terms of the "bare" bias parame- 
ters c„, we coarse-grain both tracer and matter fields on 
some scale Rl- The requirement that the final expres- 
sion for observable tracer correlations be independent of 
Rl provides us with a qiiantitativc estimate of the lim- 
its of the ansatz Eq. (1). Once an i?i-dependence is 
found, we are guided to find an additional dependence of 
the tracer density on a "regional" property of the mat- 
ter density field which absorbs the dependence on Rl 
through renormalization. We show this explicitly in two 
cases: a bias with respect to V^(5, previously found specif- 
ically for peaks of the density field, which absorbs the 
i?L-dependence induced by the smoothing of the matter 
correlation function over the coarse-graining scale; and a 
bias with respect to the variance of the small-scale den- 
sity field, which has to be introduced in the presence 
of primordial non-Gaussianity. In the latter case, the 
renormalization absorbs the CiC2(^i:,(l)(5£(2)) term into 
an overall scale-dependent bias coefficient, similarly as 
shown in [18]. However, in our approach we find a clear 
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physical interpretation of this coefficient, which preserves 
the connection to the Gaussian bias derived in [13, 14]. 

In order to achieve these results efficiently, we adopt 
several simplifications. First, we work in real space rather 
than Fourier space, since this is where the coarse-graining 
and "separation of scales" is defined. We thus do not ad- 
dress effects such as stochasticity and exclusion, which 
are restricted to very small physical separations (while 
they in general affect Fourier-space correlations at all k) . 
Furthermore, our approach is primarily intended as be- 
ing applied in Lagrangian space. However, we will not 
restrict the treatment to a linear or Gaussian density 
field. Thus, our results are applicable to biasing with re- 
spect to a nonlinearly evolved matter density field as well, 
though, in that case, further non-local biases should in 
general be considered [19-21]. We also assume that the 
tracer abundance solely depends on the total matter dis- 
tribution, rather than the baryon and cold dark matter 
density separately; the approach can easily be general- 
ized to deal with the two-component fiuid. Finally, we 
restrict ourselves to two-point correlations of tracers. 

Implicitly, we will work in synchronous- comoving 
gauge throughout. That is, all comoving observers on 
a constant-t slice share the same proper time (at linear 
order), and thus see a Universe of equal age. For a dif- 
ferent choice of time slicing, tracer correlations in gen- 
eral receive contributions from the different evolutionary 
stage of different regions as well [22] (see also [23-26]). 

The outline of the paper is as follows. In Sec. II, we 
describe the basic approach assuming purely local bias- 
ing as in Eq. (1) (although wc never actually use this 
relation) and Gaussian initial conditions, and derive the 
expression for tracer correlations in terms of PBS bias 
parameters. In Sec. Ill, we show how a bias with respect 
to V^iJ naturally appears in this approach. Finally, in 
Sec. IV we consider the case of non-Gaussian initial con- 
ditions. We conclude in Sec. V. The appendix contains 
the derivation of the key result of Sec. II for general non- 
Gaussian initial conditions, the extension of the treat- 
ment in Sec. Ill to higher order, and detailed derivations 
of some relations used in the main text. 



The small-scale density field, which we will consider in 
Sec. IV, is then defined as Ssi^) = (5(x)-(5l(x) [Eq. (86)]. 
The number of tracers (orange dots in Fig. 1) within this 
region U is simply given by the weighted sum, i.e., the 
discretized analog of Eq. (4), 

n^(x) = ^I^i(xi-x), (5) 

i 

where the sum runs over all tracers in the (idc^alizcd) 
sample, and Xj is the position of tracer i. The expecta- 
tion value (n/i), estimated by averaging over N different 
regions U and letting — > oo, is equal to the cosmic 
mean of the abundance of tracers. It can be measured 
(with some uncertainty) either observationally or in sim- 
ulations for any given tracer. 

We can now implicitly define a function Fft„i(<5L;x) 
through 

n?»(x) =F/.,l(<5l(x);x). (6) 

The dependence of Fh^^ on x denotes any departure, or 
"scatter" , of the tracer number density from a determin- 
istic relation hh{x) = hh[SL{^)]', by definition, this scat- 
ter is equivalent to the dependence of fih on the small- 
scale fluctuations Sg in the given region. The key assump- 
tion we will make below is that the correlation of this 
scatter with large-scale perturbations (in particular on 
the scales we are measuring correlation functions) is neg- 
ligible. Then, the scatter will add noise to the measure- 
ment, but will not contribute to the expectation value of 
correlation functions on large scales [27]. If this assump- 
tion breaks down, it is straightforward to include other 
properties of the density field as arguments of which 
will be the subject of the following sections. The PBS ar- 
gument we will apply below will allow us to derive the 
statistics of the tracer without any explicit knowledge 
of the function Fh,L- As indicated by the notation, the 
function Fh^L will depend on Rl. 

We can formally expand Eq. (6) in a Taylor series, 



II. PEAK-BACKGROUND SPLIT AND 
TRACER CORRELATIONS 

Consider a filter function Wi,(x) of characteristic size 
i?L, normalized to unity in 3D space. In the fol- 
lowing we assume that the filter function is isotropic, 
Wl = Wl{\x\). This is a natural assumption since any 
anisotropy would correspond to introducing preferred di- 
rections. We define the filtered (coarse-grained) density 
field in terms of the full density field ^(x) through 

5l(x) = J d^y Wl{x - y)d{y). (4) 

We can think of 5l(x) as the average density within a 
region U of size Rl centered on x (circles in Fig. 1). 



nh{^) = T.^,Pt:liO;^){SL{^)r^ (7) 

n=0 

where fI^^^{0:x) denotes the n-dcrivativc of Ffi i{SL,x.) 
with respect to 6l evaluated at position x and at 6^ — 0. 
We now take the expectation value of Eq. (7) in order 
to obtain an expression for the mean number density of 
tracers, our first observable. Our assumption of negligible 
correlation between and the scatter encoded in the 
explicit x-dependence of Ff^^L implies that the two factors 
in each term of Eq. (7) are independent random variables 
(cf. the Poisson clustering model [4]): 

(f(2(0;x) [6U^)r) = (f(;2(0;x)) ([^..(x)]") (8) 



We then obtain 
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where we have defined 



{Fh,Lm 



(9) 
(10) 
(11) 

(12) 



dropping the argument x since the ensemble average is 
independent of location due to homogeneity. The c„ de- 
pend on but for clarity we will not indicate this de- 
pendence explicitly. Further, as in Sec. I, 
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\WL{k)\^P{k). 



(13) 



Note that by definition, {5l) = 0, and cq = 1. In 
the limit that i?^ — ?> oo so that ctl — >■ 0, we see that 
{nil} = {Fh.Li^))i the expectation value of the function 
Fh,L at the background density. For finite values of Rl 
however, (hh) receives contributions from the variance 
cr| and higher order moments of the density field coarse- 
grained with Wl- This just says that for finite regions 
{Fh.Li^)) does not give the cosmic mean of the tracer 
abundance, (n^). This is commonly phrased as a non- 
zero zeroth order bias parameter given by (fi^) / (-F)i,l(0)) 
so that the cosmic mean is recovered upon taking the en- 
semble average. 



A. Correlations 

We now turn to the correlation function of tracers. 
If we measure the correlation function at separation r, 
we clearly need Rl < r in order to avoid large effects 
of the coarse graining. However, as discussed in Sec. I, 
the precise value of the coarse graining scale should not 
have an effect on the final expression for the correlation 
function. We will deal with the effects of the coarse- 
graining on the tracer correlation function in Sec. III. 

In terms of the coarse-grained densities fih, the sim- 
plest estimator for can be written as 



E 



^h(Xi)nft(Xj) 



- 1. 



(14) 



r<|xi — Xj |<r+Ar 



Here, is the mean observed abundance of tracers, 
which can be defined as fifi — hfi(x.i), where the 

sum runs over a large number N of random locations. We 
now let Ar go to zero, and take the expectation value of 
Eq. (14): 



iUr)) 



(n^(l)nfe(2)) 



1, 



(15) 



fih{xi) = Fh[SL{xi);xi] 




Ss{x) = d{x) - 6l{x) 




FIG. 1: Sketch of the separation of the density field (blue, 
thin line) into large-scale part 5l (red, thick line) and a small 
scale part Sa (Eq. (86) in Sec. IV; thin black line below), 
via an arbitrary coarse- graining scale Rl- The tracer density 
coarse-grained on scale Rl (circles) is described by the func- 
tion _Fh,L(5_t;x) [Eq. (6)], where the explicit dependence on x 
encodes the scatter around the mean relation with Sl, which 
is assumed to be uncorrelated with Sl- 



where '1' and '2' stand for two arbitrary locations sepa- 
rated by a distance r. Following the reasoning above, the 
derivatives of Fh,L with respect to Sl at the two locations 
separated by r are independent random variables. 



(^^7(0;xi)F^T(0;x.)) 



(m) 



(FtiiO;^,)) (Fi5(0;x,)). 

(16) 

Using Eq. (7), we can then write Eq. (15) in terms of the 
statistics of the coarse-grained density field Sl, and the 
coefficients c„ : 
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CXD 

n.m- 



-{Slil)Sf{2)) 



n,m— 



n!m! 



1. 



(17) 



Similarly, we can define the tracer-matter cross- 
correlation function (assuming that Sl is observable 
somehow) and obtain its expectation value: 



6im(f) 



{ihm{r)) = 



E 



nh(xi)SL{:x:j) 



r<|x^— Xj|<r+Ar 

(%(1)'5l(2)) 



nh 



(18) 



^oo c. 



Er=i|('^E(i)'^L(2)) ^ 



(19) 



These expressions involve sums over moments of the den- 
sity field, which contain disconnected pieces such as {S'j), 
multiplied by the bare bias parameters c„. The conver- 
gence of this sum of i?L-dependent coefficients multiply- 
ing i?i-dependent disconnected moments is apparently 
only assured in the general case if (t£, and higher mo- 
ments of the density field, are much less than one, since 
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terms of order appear at arbitrarily high n. On 
the other hand, a physically reasonable perturbative bias 
model for ^/i(r), ^hmif) should converge as long as the 
connected matter correlators are much less than one, in- 
dependently of the choice of the fictitious coarse-graining 
scale. Thus, our goal is to reorder the sum in Eq. (17) 
and Eq. (19) into a sum of iJ^-independent coefficients 
6]v multiplying only powers of connected matter correla- 
tors. This can be seen as a rcnormalization of the "bare" 
coefficients c„ into "renormalized" bias parameters 6jv. 
We will see below that they have clear physical signifi- 
cance. 

We can expand the correlators into connected parts 
(cumulants) as follows. The correlator appearing in the 
cross-correlation Eq. (19) can be written for any statisti- 
cal field as 

(5£(1)542)) = (n) {Sr''HS^aW))c, (20) 
jv=o ^ ^ 

where the subscript c denotes connected correlators. 
Thus, Eq. (19) can also be written as 

^ oo / \ 

n=l ■ Ar=0 ^ ^ 

(21) 

where we have defined 

oo 

^^^^(<5E). (22) 

71=0 

The density field correlator in the auto-correlation 
Eq. (17) is more complicated. As shown in App. A, 

(SiST) = EE (fc) {Si) h) {sl){sr'sT-%.i , 

fe=o ;=o ^ ^ ^ ^ 

(23) 

where (•)nzi denotes a disconnected correlator that, when 
expanded into cumulants, does not contain any zero-lag 
pieces, i.e. no factors that asymptote to a constant as 
r ^ oo (see Eq. (A12) in App. A for a mathematical 
expression of this definition). We will restrict to the 
Gaussian case here, and postpone the discussion of non- 
Gaussian density fields to Sec. IV and App. A. The no- 
zero-lag requirement entails that (5l(1) and (5l(2) have to 
appear in equal powers in "nzl" correlators. Defining 

a(r) EE {5UI)Sl{2)), (24) 

we obtain 

(^£(l)5r(2))nzi = n! [^(r)]" Snm- (25) 

The factor n\ comes about because there arc n! ways to 
contract 2n factors of 5l {n of (5^,(1) and Sl{2) each) into 
a product of n correlation functions. Further, 

{6l) = {n-l)\\al (26) 



for n even, and zero for n odd. This yields 

min(n,m) 

(52(l)<5r(2)> °'"'='^'' E {n-k-l)U{m-k-l)\\ 

k: n—k^ m—k even 

xQ(7)fc!ar™-^'^[a(r)]'=. 

(27) 



B. PBS and bias parameters 

So far, all we have done is divide the universe into ficti- 
tious regions, and describe the number density of tracers 
in regions in terms of a function Fh,L [Eq. (6)]. We were 
then able to formally express the correlations of tracers 
in terms of the statistics of the matter density and the 
derivatives of the function Fh^L, all of which depend on 
the coarse-graining scale Rl. 

We now turn to the peak-background split argument, 
and the definition of the PBS bias parameters with re- 
spect to density. The argument can be summarized as fol- 
lows: if the description of the clustering of tracers solely 
through their dependence on Sl is sufficient, then the ex- 
pected abundance of tracers in a region U characterized 
by a coarse-grained overdensity Sl = D is sufficiently 
well approximated by the average abundance of tracers 
(hh) in a fictitious Universe with modified background 
density 

p' = p(l + D) , (28) 

where p is the actual background density. 

The advantage of this approach is that we only need 
a prediction for (n/i) as function of the background den- 
sity p' to calculate the statistics of tracers; no knowledge 
of the function Fh,L is necessary. Note also that this is 
directly connected to the derivation of bias (linear bias 
in that case) in the relativistic context presented in [22]. 
Specifically, we are working in the synchronous gauge 
where all space-time points on an equal-coordinate-time 
hypersurface share the same cosmic age. Correspond- 
ingly, when calculating {fih) for varying p' it is crucial to 
keep the age of the Universe fixed. 

Thus, we now consider the case where we perturb the 
background density by Ap ~ Dp, where D is an infinites- 
imal parameter. Thus, in a region with overdensity 6l 
the matter density is perturbed to 

PL = p(l + Sl) ^ p(l + Sl) + Ap = p(l + Sl+D) (29) 

Note that we add a fixed amount of uniform matter den- 
sity everywhere; we do not rescale the local matter den- 
sity yO by 1 + -D, which would also amplify the fluctuations 
S. We can obtain the average number density of tracers 
(hh) (more precisely, the expectation value of the esti- 
mated mean number density in some volume) in such a 
Universe from the expansion in terms of coarse-grained 
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Sl, Eq. (11): 



{nh)\D 



(F,,L(0)>X:^(('5i+^r>, (30) 

71=0 



where Fh,L and c„ both refer to the Universe with back- 
ground density p, i.e. D = 0. 

Let us now define the peak-background split bias pa- 
rameters 6jv {N > 1): 



6jv 



1 



{nh)\D= 



D=0 



Using Eq. (28), we can also write this as 



6jv = 



{nh) dp 



N 



(31) 



(32) 



where the derivatives are evaluated at the fiducial value 
of ~p. It is worth emphasizing the difFcrcncc between those 
bias parameters and the c„ defined in the last section: the 
6]v quantify the response of the cosmic mean abundance 
of tracers to a change in the background density of the 
Universe; specifically, they do not make any reference to 
the regions U, or the scale Rl- The c„ on the other hand 
quantify the average response of the abundance of tracers 
within a region hi to a. change in the average density 
5l within that region, evaluated at 5l = 0; they thus 
necessarily depend on the filter function Wl and scale 
Rl- Further, there is no "scatter" in the defining relation 
Eq. (31) for the PBS biases, although the values for bN 
measured in reality will clearly have a finite error bar 
as we can only approximate this relation within a finite 
volume. 

The 6jv are closely related to the resummed bias prop- 
agators defined in [11] [see Eqs. (83)-(84) there], while 
the bare bias parameters c„ correspond to the bare prop- 
agators [Eqs. (l)-(2) in that paper]. 

Using Eq. (30), we can derive an algebraic relation 
between 6jv and c„: 



bn = 



1 



N ^ n\ {n-N)\ 



/ cn—N\ 



(33) 



By reordering the sum in Eq. (21), we immediately see 
that the tracer-matter cross-correlation in terms of the 
PBS bias parameters is given by 



JV=1 



(34) 



This is the well-known bias expansion of the tracer- 
matter cross-correlation function. However, note that 

the matter correlators appearing here arc the connected 
correlators. In particular, in the case of a Gaussian den- 
sity field we obtain 



(35) 



i.e. the tracer-matter cross-correlation function is sim- 
ply proportional to the linear matter correlation func- 
tion (the same conclusion was reached by [28], who 
only considered Gaussian density fields). Similarly, us- 
ing Eq. (33) and Eq. (27) we can re-express the auto- 
correlation [Eq. (17)] as 



iUr)) 



OO 



bjsf bM 



(5f(l),5f (2)) 



M/ 



^ N\M\ 



nzl • 



(36) 



Even though Eq. (27) assumes a Gaussian density field, 
Eq. (36) is in fact valid for a general non-Gaussian density 
field. The proof in this case requires somewhat more 
effort and is given in App. A. Note that there is no 60; 
the expressions Eqs. (34)-(36) only involve terms with 
6jv for AT > 1. In the Gaussian case, Eq. (36) further 
simplifies to 



m 

JV=1 



(37) 



Eqs. (34) (37) achieve the desired result: an expan- 
sion of the tracer correlation function in terms of R^- 
independent bias parameters which multiply powers of 
the matter correlation function ^Lif) (or, more generally, 
no-zero-lag correlators). Thus, this expansion converges 
if and only if the finite-lag matter correlators are small on 
the scale r, which is what we expect from a physical bias 
expansion. On the other hand, in the bare bias expan- 
sion, Eq. (17), terms of order cr£ appear at every succes- 
sive higher order, suggesting that we need to coarse-grain 
the density field on quasi-linear scales in order for the ex- 
pansion to be perturbativcly valid. This of course would 
be disastrous for any sharp features in the correlation 
function; for example, choosing a coarse-graining scale 
of Rl = 30 — 50Mpc//i would erase the baryon acous- 
tic oscillation (BAO) feature at r = 150Mpc//i. In the 
PBS bias expansion, there is no need to choose a coarse- 
graining scale this large. Rather, the validity of the the 
result Eqs. (34)-(37) is determined by the requirement 
that any residual dependence on JS^ be negligible. 

Another important property of this expansion is that 
the same PBS bias parameters describe both the matter- 
tracer cross-correlation and the tracer auto-correlation, 
which is what we expect from a physical bias expansion 
(the corresponding statement in Fourier space is compli- 
cated by the small-scale effects such as stochasticity, shot 
noise and halo exclusion which contribute to the tracer 
power spectrum at all k). 

One crucial advantage of this approach is that we 
have an indicator for when the underlying assumptions 
break down: if evaluation of Eq. (36) shows that the 
result is in fact i?z,-dependent, then we know that the 
underlying assumptions, in particular the description 
n;i(x) = Fh^L{SL{^)), break down. One then has to iden- 
tify the physical reason for this i?L-dependence, and add 
a dependence of the tracer number density on additional 
parameters which will absorb ( "renormalize" ) the Rl- 
dependence. We will see two important examples of this 
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in Sec. Ill and Sec. IV. First however, wc will derive 
the PBS biases for the widely considered case of univer- 
sal mass functions, and then illustrate the approach on a 
concrete example of biased tracer: regions above thresh- 
old. 



C. PBS biases for universal mass functions 

The mean abundance of tracers such as dark matter 
halos of some mass M* is often parametrized in the form 

nh=P f{vc) J* (38) 
5c In cr« 



J* 



din i?. 



(39) 



where a, is the variance of the linear matter density 
field on scale R^, is related to the mass M» through 
= 47r/3pi?i^, and 8c is the linearly extrapolated 
threshold for collapse. Further, /(i^c) is in general an 
arbitrary function of v^- The Jacobian J* is present in 
order to convert from an interval in cr* to a mass interval. 
Eq. (38) is referred to as "universal mass function" and 
was originally motivated by the excursion set formalism 
[29] . It is a special case of a more general description of 
mean tracer abundance we will consider in Sec. IV E. 

In order to derive the bias parameters Eq. (31), we 
need to know how hh changes under a change in the 
background density of the Universe [Eq. (28)]. Since we 
work in the Lagrangian picture, we will ignore the triv- 
ial dependence through the p prefactor in Eq. (38). The 
threshold is defined as the fractional overdensity a re- 
gion miist have to collapse^ at a fixed proper time to- In 
an Einstein-de Sitter Universe, a spherical perturbation 
with a mean initial fractional overdensity 6c ~ 1.686, i.e. 
with p{< R, t) = [1 + a{t)Sc] pit) average interior density 
for a{t) <§; 1, collapses at o = 1. The same reasoning 
also holds for more general expansion histories, where 
5c assumes other values. Since the evolution of such a 
perturbation is independent of the external Universe (by 
BirkhofF's theorem), a perturbation of the same physical 
density pc will collapse at the same proper time in a Uni- 
verse with perturbed background density p' =~p{l + D) 
as well. The significance = 5c/ cr* = {pc — p) /5prms 
quantifies how rare fluctuations above a physical density 
threshold pc = {1 + 5c)p are given the RMS fluctuation 
amplitude ^prms = cr^p. Clearly, if we add a uniform 
matter density component Dp, the critical overdensity 
changes to 



Pc-p' = {l + 6c)p-{l + D)p={6c 

Thus, the significance is modified to 

^/ ^ Pc-p' ^ 5c- D 
" a^p cr* 



D)p. (40) 



(41) 



^ Since General Relativity is scale-free, this threshold is indepen- 
dent of the size and enclosed mass of the perturbation. 



For a mass function of the form Eq. (38), changing the 
background density is thus equivalent to changing 5c — >■ 
5c — D. Eq. (38) and Eq. (31) thus immediately yield 



6jv = 



(-i)^a^(%) (-1)^ 1 d^f{uc) 



{hh) 85^ 



(42) 



This is the widely known expression for the peak- 
background split bias parameters, which in our approach 
is a special case of Eq. (31). 



D. Application to regions above threshold 

We now turn to a simple example of tracer for which an 
exact expression of the tracer correlations is known. Pre- 
cisely, we define our tracer to be a region where the den- 
sity field 5 is above a fixed threshold 5c = t'cf* , where cr* 
is the RMS fiuctuation of the density field and Vc an arbi- 
trary fixed parameter. The density field can be thought 
of as smoothed on some scale r*; however, since for our 
purposes this scale is irrelevant, we will not make this 
smoothing explicit in our notation in order to avoid con- 
fusion. Note also that we do not make any assumption 
about fc, such as the "high peak limit". In the present 
case, unlike peaks of the density field, the tracer popula- 
tion is not a (countable) point set. Rather, tracers cover 
a finite volume, and the tracer number can be defined as 
a continuous field 



iVft(x) = e(5(x) - VcT^). 



(43) 



The number density of tracers as defined earlier in this 
section is then given by 

n^x) = j dV^/.(y)W^L(x - y). (44) 

We first review the exact approach to the clustering of 
such "tracers", and then investigate the PBS prediction. 
Note that essentially all these results have already been 
derived in [11, 30-35]. However, we review it here in light 
of the discussion presented above. 



1. Exact calculation for a Gaussian density field 

The mean "number density" {hh) of the tracers defined 
above is simply given by the fraction of the total volume 
that is above the threshold Vcf^* , 



(45) 



where Piivc) is the probability that the density field at 
a random location is larger than fcC*- If the underlying 
density field follows Gaussian statistics with variance cr, , 
Pi is given by 



1 



'da;e--V2= lerfc(^-^ 



(46) 
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The exact expression for the two-point correlation func- 
tion of our "tracers", ^^(r), is then given by the proba- 
bility of finding two peaks P2 separated by r, relative to 
the random probability (Sec. Ill B of [34]) through 



P2{iyc;r) 

erfc 



- 1 

V2) 



1 N 



(47) 



Here, ^(r) is the two-point correlation function of the 
underlying density field smoothed on the scale i?*, so 
that ^(0) = al 



2. Peak-background split calculation 

Eq. (36) gives the tracer correlation function in terms 
of the PBS bias parameters and powers of ^i,(r). Note 
however that given the explicit relation between hh and 
6 through Eq. (43) in this simple example, the; division 
into regions U is a, purely conceptual device here, and 
we can always set Rl to be equal to or smaller than the 
smoothing scale adopted in the thresholding approach. 
Hence, we will drop the subscript Rl below. 

The PBS bias parameters are defined by Eq. (31). 
Given our definition of tracers as regions where S{x) > Sc, 
or equivalently, p(x) > p(l -|- 6c), we see that a fractional 
change D in the background density is equivalent to a 
change in the threshold 6c- 

p^p{l + D) ^ 6c^6c-D. (48) 

Hence, the PBS bias parameters are given by 

^ (-l)^9^(n,) ^ (-l)^ 1 d^P, 

{hh) 86^ <7f Pi{uc) dv^ ■ ^^^^ 

This of course can also be derived by noting that the 
"abundance" of regions above treshold Eq. (45) is a spe- 
cial case of universal mass functions (Eq. (38) without 
the Jacobian factor which is irrelevant for the ^at), so 
that Eq. (42) applies. By using the generating function 
of the (probabilists') Hermite polynomial 

H^{x) = (-l)"e^^/^|; (e-V2) (50) 

we calculate the n-th derivative of Pi (n > 1) as 



= (-l)"if„_l(l.e) 




(51) 



Thus, wc can explicitly write the PBS bias parameters 
for our tracers: 



3JV 



erfc I 

\V2 



-HN-iiuc). (52) 



Inserting this into Eq. (37), we immediately obtain the 
PBS prediction for the correlation of thresholded regions: 



a(r-) 



'''' ' 72) 



-2 
e~ 



]V=1 ■ * 



(53) 



We see that this agrees with the direct calculation, 

Eq. (47). A mathematically similar derivation was pre- 
sented in [34]. The difference is that in [34], we used 
Eq. (49) and Eq. (47) to infer the general relation 
Eq. (37). Here, wc are simply illustrating how the inde- 
pendently derived Eq. (37) applies to the case of thresh- 
olded regions, a case where we know explicitly the func- 
tion P;i,l((5l). As proven in the previous section, Eq. (37) 
and the much more general Eq. (36) apply to any tracer 
as long as the dependence of the tracer density on other 
quantities apart from the matter density can be ne- 
glected. 

One alternative to the simple local bias expansion 
Eq. (1) in the context of thresholded regions is to ex- 
pand Eq. (43) in terms of Hermite polynomials, as done 
in [11, 33, 35] (see also [2, 4]): 



6(1/ - uc) 



n=0 



an(yc)iin(y) 



(54) 



where i/(x) = (5(x)/(t*, and 



{^.^e-^'J^Hr^.M n>l 



(55) 



The bias parameters they obtain are exactly equal to 
our renormalized PBS bias parameters, Eq. (52). The 
reason for this is that, in the Gaussian case, the Hermite 
expansion ensures that no disconnected pieces remain in 
the correlation function expression Eq. (17). Specifically, 
denoting = i/(xj), we have 



{Hr,{u,)Hm{iy2)) = 



1 

2^ 

X exp 



V 0-2 dPidu2 



dV2 Hmiv^) 



1 °° 1 

— V — 

27r ^ m 

JV=0 



N 



J du^Hn{vi)HN{vi)e-''"'^ 



X / dv2H^{u2)HN{y2)e"'i''^ 

e(ri2) 



(56) 
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The Hcrmitc expansion is thus an elegant way of directly 
obtaining rcnormalized bias parameters in the case of 
thresholding in a Gaussian density field. However, the 
additional contributions obtained in the non-Gaussian 
case spoil this renormalization of all zero-lag terms, as 
we will see in Sec. IV. 



III. SMOOTHED CORRELATION FUNCTION 
AND CURVATURE BIAS 

Above wc explained that the expression of the tracer 
correlation in terms of PBS bias parameters and con- 
nected matter correlators [Eq. (36) and Eq. (37) for the 
general and Gaussian case, respectively] should be nu- 
merically insensitive to the value of the; coarse-graining 
scale R]^. Further, a significant i?i-dcpendcncc indicates 
the break-down of our assumption that the tracer density 
is a function only of the local coarse-grained density. 

For r ^ Rl, and as long as ^(r) is smooth (e.g., a 
power law), the smoothed version of ^(r), ^^(r), will not 
differ significantly from ^(r). However, if ^(r) has some 
features on a scale Sr <C r, such as the BAO feature with 
dr ~ 20Mpc/ft,, then the condition for i?i-independencc 
becomes much more restrictive: Rl ^ Sr. We now show 
how the i?i,-dependence induced through ^z,(r) can be 
cured. 

For a general isotropic filter function (see Eq. (Bl) in 
App. B) the effect of smoothing on the correlation func- 
tion ^(r) can be perturbatively described through 



a(r) 



(2^ 



\WL{k)fP{k)e 



ikr 



: JA±^l_2Rie + 0{k'))P{k)e 
--ar) + 2RlV^ar) + 0(W^C{r)), 



ik.r 



(57) 



by suitable definition of the parameter Rl. In App. B, 
we give the general expansion of WL{k) [Eq. (B2)] and 
the expansion of ^L{f) in terms of derivatives of ^(r) 
[Eq. (B3)]. Thus, if ii|,V^^ is comparable to ^, our 
requirement of i?L-independence does not hold. In our 
approach, this signals a breakdown of the underlying as- 
sumption that tracer statistics can be described purely by 
their dependence on the local matter density. Instead, let 
us assume that the local number density also depends on 
the coarse-grained Laplacian of the density field, i.e. the 
curvature: 



n^x) = F^,i(,5i(x); V'5l(x);x). 



(58) 



The Laplacian is the lowest order term in derivatives of 

Sl, because a dependence on the gradient of Sl would 
imply a preferred direction'^. In general, we now have 



^ Terms such as (VSl)^ could however appear at second and 
higher order. 



to perform a bivariate expansion of the function F^.l in 
6l and V^Jl- Let us for now restrict to lowest order 
to keep the treatment clear, and consider the Gaussian 
case. The expansion to higher orders in derivatives of 
5l is described in App. B. The tracer auto-correlation 
becomes 

a(r) = cI{Sl{1)6l{2)) + 2ciCv25(5i(l)V2(5i(2)) 
+ O(V^C) 

= cl [^(r) + 2i?iv2^(r)] + 2c^c^2sV^^{r). (59) 
Here, we have defined 



Cv2 5 



1 



OF, 



h,L 



5l=0,V^5l=0 



(60) 



Eq. (59) is again phrased in terms of (in general) discon- 
nected matter correlators and i?L-dependent bare bias 
parameters. We now need to introduce a iJ^-independent 
PBS bias parameter for V^(5 as defined in Sec. II B for 
the density itself. We would like a transformation where 
the Laplacian of the density perturbation shifts by a con- 
stant: 



V2(5a(x) = V2<5(x) + 



a 

£2' 



(61) 



where a is a dimensionless small parameter, and we have 
added a length scale £. This corresponds to 

5„(x)=5(x) + ^(x2 + A.x + C), (62) 

where A and C are constants and the center of the region 
considered is chosen as the origin. We are not interested 
in adding a gradient to the density field and hence set 
A = 0. Note that Eq. (62) is only defined for a region 
of finite size (e.g., a simulation box), so that (n/i) in the 
following is to be considered as an ensemble average over 
many such finite regions. We will also set C = so that 
(5(0), at the center of the region considered, is unchanged 
(a constant shift in S such as described by C just corre- 
sponds to the density bias transformation of Sec. II B). 
Thus, 



a 



^„(x) = ^(x)-F^x-. 



(63) 



We can now defined a (renormalized) PBS bias parameter 
through 



(n/i) da 



a=0 



(64) 



Defined in this way, the scale (. will disappear o\it of the 
final expression for the tracer correlation function (note 
that 6v2(5 has dimension length squared). In order to de- 
rive the relation between 6v2^- and the c„, we need the 
behavior of both (5l(x) and V^(5l(x) under the transfor- 
mation Eq. (63). The latter is immediately obtained as 



V25L,a(x)=V25L(x) 



(65) 
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We are interested in the change of 5l near the origin 
(the center of the region U). In analogy with Eq. (57), 
we obtain 



'JL,a 



(0) = J d'yWLiy) 



(66) 



Thus, Eq. (58) yields at lowest order 

&V2^ = 



ciRj 



da 



da 



(67) 



Now, we can write the generalization of our previous re- 
sult, Eq. (37), at lowest order, assuming we absorb the 
smoothing effect on ^(r), as 



a(r-) = ^?C(r-)+2&iM5V'e(r-), 



(68) 



where the factor of 2 comes from the two permutations 
when writing down all mixed no-zero-lag correlators be- 
tween 5 and V^(5. Using Eq. (67), this is equal to 

Ur) = c^ar) + 2 [ciRl + c.c^.s] V'e(r) + 0{^^^{r)), 

exactly matching the result of Eq. (59) . Thus, by intro- 
ducing a dependence of the tracer density on the Lapla- 
cian of the density field, and a corresponding PBS bias 
parameter, we arc able to absorb the effects of the coarse- 
graining on the correlation function. Moreover, in App. B 
we show that this continues to arbitrary powers of deriva- 
tives, and that the PBS bias parameters can entirely ab- 
sorb the smoothing effect on [Eq. (B16)]. Specifically, 
up to order V''^(r), we obtain 



a(r)=6?e(r) + 26i6v25V2C(r) 



(69) 



There are two important implications of this result. 
First, our approach, which does not make any assump- 
tions on the tracers themselves, gencrically predicts the 
existence of a bias with respect to V^(5, which in A;-space 
corresponds to a scale-dependent bias oc fc^. One can in- 
terpret this as the statement that the tracer density is in 
general not a truly local function of the matter density, 
but depends on the matter distribution within a finite 
region whose characteristic scale is given by V^'v^a- 

Second, the definition of the PBS bias parameter fov^s 
[Eq. (64)] has a clear physical interpretation: it corre- 
sponds to the response of the tracer number density to 
a uniform shift in the curvature of the density field (and 
shifts of higher derivatives in the general case, Eq. (B4)). 
Below we will show how this bias can be evaluated for an 
analytical example, peaks of the density field. A quanti- 
tative test on N-body simulations will be the subject of 
future work. 



Beyond linear order in the matter correlation function, 
one in principle has to expand the tracer density in a 
multivariate bias expansion of V^"(5i (n > 0). Although 
we have only shown that this expansion removes the Rl- 
dependence contained in ^L(r) at linear order in ^, we 
expcict this to be the case for higher orders as well. In 
practice, the suppression of fov^^^V^"^ compared to ^ 
ensures that one only needs to keep a finite number of 
terms. 



A. Connection with the peak model 

In the peak model, we identify large-scale structure 

tracers with discrete peaks of the density field above some 
threshold Sc- These distinct peaks constitute a point set 
(in contrast to the regions above threshold considered in 
Sec. II D). For a Gaussian density field, it is possible to 
calculate the two-point correlation of these peaks exactly 
[36, 37]. It is well known that peaks exhibit a scale- 
dependent bias in Fourier space cx k"^ which is equivalent 
to a bias with respect to V^(5 [36-38]. In this section, we 
show how the approach outlined in the previous section 
relates to this model. 

As in Sec. II D, we smooth the density field on a scale 
i?* and define the local significance v{:x) = 5(x)/(t*. 
Note that the smoothing scale here is physical, unlike 
the fictitious coarse-graining scale adopted in the renor- 
malization approach described above. Apart from i?,, 
which is usually identified with the Lagrangian radius of 
the lialos c;onsidered and is irrelevant for the discussion 
here, the peak model by itself does not involve any other 
coarse-graining scale. We begin by outlining the deriva- 
tion of the mean number density of peaks, which follows 
App. A of [39] , but will use the slightly different notation 
of [36, 37]. 

The location Xp of a peak is defined through con- 
straints on the overdensity S (or equivalently u), its gra- 
dient rj = V5, and its Hessian = didjS: 



T7(xp) =0 
Ai(x), A2(x),A3(x) >0, 



(70) 



where Vc = Sc/cr^, is the scaled threshold, and arc the 
eigenvalues of Qj- Let V denote the 10-component vector 
consisting oiv,ri, and the six independent components of 
Qj . The differential (in terms of v) mean number density 
of peaks is then given by 

npk(i'c)(x j 5D{iy - Vc)5d{ii)Q{\i)Q{X2)Q{\z) 



Q(V) = -V^M-^V, 



(71) 



where M is the covariance matrix of V. The most difia- 
cult part of the calculation is deriving the integration re- 
gion and measure of V, but we will not need to deal with 
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this explicitly as wc arc only interested in how npk trans- 
forms under a change in Inp and under the transforma- 
tion Eq. (61). Following Bardeen et al. [39] we introduce 
new variables u,y,z {u is their x), where in particular 



u = —- 



0-2 



0-2 



4 = ii^'Sf), (72) 



corresponds to minus the curvature of the density field 
scaled to unit variance. The log-likelihood Q then be- 
comes 



2Q 



2 



1-72 



+ 2Q 



where 



^v, and 7 = {vu) 



(73) 



(74) 



is a scaled spectral moment quantifying the correlation 
between v and u. 2Q^ is the log-likelihood of y, z, t], 
and the other components of Cij, and is independent of 
and u. 

Bardeen et al. [39] then find for the mean differential 
peak number density 



1 



(27r)2ii3 



Go(7,7^^c)e-''='/^ 



(75) 



where i?i is the characteristic radius of a peak [Eq. (19) 
in [37]], and 



exp 



{u-u,y 

Goi,,u.) = J^ dufiu) ^^^^^^ 



(76) 



Here, f{u) is a function encoding the integral over y and 
z which accounts for the asphericity of the peak profile 
[39]. 

Going back to Eq. (73), we see that adding a uniform 
density component corresponds to a change in the thresh- 
old as described in Sec. II D, whereas the other variables 
u, y, z are not affected. Thus, the density bias parameters 
are given by 



6jv 



(-1)^ a^np k 
fipk d{Sc) 



N' 



(77) 



where is the height of the density threshold. These 
PBS bias parameters precisely agree with the coefficients 
of ^(r) and [^(r)]^ derived from an explicit computation 
of the peak 2-point correlation function [37] . The relation 
Eq. (77) was already pointed out in [15]. 

More interesting in this context is the derivation of the 
bias with respect to V'^S. We consider the transformation 
Eq. (61) in a region of finite size £, so that the effect on the 
density and v is negligible if a <i; 1. Thus, all components 
of the vector V arc unaffected with the exception of u. 
Instead of following a Gaussian of mean zero and variance 



of 1, u is now Gaussian-distributed around (recall the 
minus sign in the definition of it) 



a 

"02^' 



(78) 



with unit variance. Thus, the log-likelihood Eq. (73) 
changes to 



2Q(a) 



1 



a 



1 - 72 V fT2^2 

Clearly, this is equivalent to changing 



M* - 



a 
02^' 



2Q+. (79) 



(80) 



The remainder of the calculation of fipk follows through 
as before, since it is independent of the distribution of 
u. Wc thus obtain for the renormalized curvature bias 
[Eq. (64)] in the peak model 



npk da 



1 a^pk 



a=0 



npkO'2 du 

5Go(7)W*) 



Go(7,7J^c)i72 



= — (GoCr2 



f° 

Jo 



du 



u — 'yu, 
1-72 



2(1-7^) 



1 

CT2 



u - 71/c 



1 



(81) 



where u is the mean peM,k curvature (i.e. the intc\s;ral of 
u times the integrand of Go, normalized by Gq and eval- 
uated at 71'c)- Eq. (81) is precisely the scale-dependent 
bias parameter 601 found in [36, 37], except for a minus 
sign which arises from the fact that these authors defined 
601 as the linear bias associated with = —"V^S. The 
biases with respect to higher powers of V^S are obtained 
by generalizing Eq. (81) to 



(82) 



In particular, we find for the PBS bias with respect 



5(V25)2 



(G 



oai)-i / 
Jo 



°°dJ"-^':'''-'!-^''/(n) 



2\2 



2(1-7'') 



VM-L - 7^) 



pk 



(I-7T 



1 



(83) 



in agreement with 602 as derived in [40]. Note that by 
construction, the peak model does not predict any bias- 
ing with respect to higher than second derivatives of the 
density field (6v2"5 = for n > 1). 
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Desjacques et al. [37] were also able to derive the scale- 
dependent bias 601 — ^b-^^s from a peak-background 
split calculation. However, they employed a conditional 
mass function, i.e. the number density of peaks given a 
(spherical) overdensity on a much larger scale Rb ^ R*, 
whereas here we derived all PBS bias parameters from the 
unconditional mass function (this is also the approach 
taken by [40]). Nevertheless, all these treatments are 
based on a similar reasoning, i.e. a long-wavelength per- 
turbation shifts the mean curvature (in the case of [37], 
the shift is correlated with the density), and lead to the 
same final result. Note that the peak-background split 
approach can be generalized to derive full expressions for 
the peak correlation functions [40], which also includes 
dependencies on more general rotational invariants, such 
as (V(5)^ and [{didj — SijW'^ /3)S]'^ . These are present in 
the clustering of peaks, even though it is not necessary 
to introduce them in order to cure the iii-dependence of 
induced by smoothing, as we have seen here and in 
App. B. 

It is important to note that the peak model is fun- 
damentally different from local bias expansions in the 
sense that it only involves a physical smoothing scale R^ , 
and no further coarse-graining. Nevertheless, the term 
a boiV^S in is equivalent to the generic b^2g term in 
^fi- Namely, 6oiV^(5 restores the contrast of the baryon 
acoustic oscillation, otherwise smeared out by the filter- 
ing in 6fo?* [36]- Moreover, in general feoi can be greater 
than ~ R^ so that the contrast of the BAO in the peak 
correlation ^pkl?") can even be enhanced relative to that 
in the unsmoothed mass correlation ^(r) (sec Fig. 5 in 
[36]). However, subsequent gravitational evolution sup- 
presses most of this scale-dependence (expected to be at 
the few percent level at the time of collapse [37]). We 
expect this to be a general feature of a V^(5 bias specified 
in Lagrangian space. 

IV. NON-GAUSSIAN CASE 

We now return to the case of a tracer which can be 
sufficiently well described by density bias, i.e. we ne- 
glect the curvature bias corrections, but consider the case 
of initial conditions that are non-Gaussian. Specifically, 
we will first focus on the case of local primordial non- 
Gaussianity, which has been shown to lead to a large 
modification of clustering on large scales; general shapes 
of non-Gaussianity will be considered in Sec. IV D. 

The derivation in Sec. II can be generalized to a general 
non-Gaussian density field, in which Eq. (36) formally re- 
tains its validity (see App. A). However, one can easily 
show that in general the resulting tracer correlation func- 
tion depends on the coarse-graining scale Rl- At lowest 
order the tracer auto-correlation becomes 

a(r) = blUr) + bMSL{l)6l{2)) + 0{St), (84) 

the second term being the leading non-Gaussian correc- 
tion. For primordial non-Gaussianity of the local type. 



and in the limit r ^ Rl, the second correlator is given 
by (see Eq. (97)) 

{Sl{1)51{2)) = 4/NL<7iC^,(r), (85) 

where is the cross-correlation between the matter and 

primordial Bardcen potential (j). Given the appearance 
of cr£, Eq. (84) is strongly i?i-dependent. This indicates 
that the description of the tracer density as a function of 
the matter density Sl alone is insufficient even on large 
scales in the non-Gaussian case. 

Instead, we need to include a dependence of the tracer 
density on the amplitude of small-scale fluctuations. This 
dependence is present regardless of the nature of the ini- 
tial conditions; however, only in the non-Gaussian case 
are there large-scale modulations of the small-scale fluc- 
tuations, due to mode coupling, whereas in the Gaus- 
sian case we were able to neglect the small-scale fluc- 
tuations in the large-scale description. In general, one 
would imagine that the abundance of tracers depends on 
the amplitude of small-scale fluctuations on a range of 
scales. However, for simplicity we will parametrize the 
dependence through the variance of the density field on 
a single scale i?*. In the local model, which we will fo- 
cus on here, this is sufficient in any case as all small-scale 
fiuctuations are rescaled equally (in the large-scale limit), 
so that the value of R* becomes irrelevant for the final 
result. 

While we focus on primordial non-Gaussianity of the 

local type here, the extension to other types of non- 
Gaussianity is straightforward (see Sec. IV D). Further- 
more, we only rely on the description of the density field 
in terms of Af-point functions, with the 3-point function 
being the lowest order non-Gaussian contribution which 
we fociis on here. That is, wc do not rely on a ficti- 
tious Gaussian field from which the non-Gaussian field is 
constructed. This is different than the approach taken in 
[16, 41, 42], where the separation of scales is typically ap- 
plied in the fictitious Gaussian field, and an application 
in the physical non-Gaussian potential is not straightfor- 
ward to implement [16]. 

We first define the small-scale density field as the local 
fluctuations around the coarse-grained fleld Sl- 

Ss{x)=54^)-Sl{x) (86) 
= J d'y[W4^ - y) _ M/^(x - y)]<5(y) 




VK,(fc)5(k)e^'^--, 



W,{k)=W,{k)-WL{k). (87) 

In Fig. 1, 5s is illustrated by the thin black line in the 
lower part of the figure. Note that as A: — >■ 0, Ws{k) oc k"^, 
i.e. the long- wavelength modes are filtered out as desired. 

This implies that the cross-correlation between Sg (x) and 
the density field (5ii(x) smoothed on some scale R goes 
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to zero as R becomes much larger than Rl: 



((5,(x)(5fi(x)) = 



j^W,{k)Wn{k)P{k) 



(88) 



and similarly for {Ss{xi)6l{x2)) if |xi — X2I '> Rl- We 
further use the notation 



a: 



2 _ 



= (Si) = 



(27r)3 



(89) 



We quantify the dependence of the tracer abundance on 
the amplitude of small-scale fluctuations through 



(90) 



where the subscript * refers to the smoothing scale 
(y*) = 0, and the factor of 1/2 is included to obtain 
expressions which conform to standard convention later 
on. In the Gaussian case, S.si'i') ^ for r » Rl, so that 
the small-scale density field and in particular have no 
large-scale correlations. In the non-Gaussian case how- 
ever, y, is in general correlated with long-wavelength per- 
turbations. Note that ((5s(1)(5l(2)) vanishes by construc- 
tion on large scales [Eq. (88)], so that it is natural to 
start the expansion with the leading term 5^. 

We now generalize Eq. (6) to explicitly include the de- 
pendence on y*, 



nft(x) = Fft,i,((5z,(x),j/*(x);x). 



(91) 



Although our approach here is formally similar to the 
bivariate local expansion in 5l and adopted in [18, 26, 
41], there is somewhat of a conceptual difference in that 
we expand fih purely in terms of properties of the matter 
distribution. The effect of non-Gaussianity, and the fact 
that it derives from a potential 0, only enter through the 
expressions for the correlators between 5l and y* here. 
The nature of non-Gaussianity thus decouples from the 
description of the tracers (which only know about the 
matter density field) in this approach. 

We can now repeat the derivation of Sec. II, includ- 
ing this additional dependence. All arguments about the 
residual scatter from the deterministic relation n/((x) = 
^?i[^l(x), y*(x)] and its negligible correlation with long- 
wavelength perturbations made in Sec. II also apply here. 
In fact, the dependence of hh{'x.) on y*(x) is a source 
of uncorrelated scatter in the Gaussian case which be- 
comes correlated with long-wavelength perturbations in 
the non-Gaussian case. This is another way of seeing 
why wc! need to introduce the dependence on explic- 
itly when dealing with large-scale non-Gaussianity. Tak- 
ing the expectation value of Eq. (91), we obtain 



^ n\m\ 



(92) 



where we have defined bivariate "bare" bias parameters 
through 



1 j d^+^Fh.L 



(93) 



We then need expressions for the various cross- 
correlations of 5l and y*. In the following, we will re- 
strict ourselves to the leading order terms, as the general 
expansion becomes lengthy. 



A. Primordial non-Gaussianity of the local type 

We will consider a density field derived from a Bardccn 
potential with non-Gaussianity of the local type. We 
will restrict our treatment to leading order in the non- 
linearity parameter /nl • At this order, the only relevant 
A^-point function is the bispectrum, 

B(ki, k2, ks) = Af(fci)X(fc2)7W(/C3)S4ki, k2, kg) 

B^(ki,k2,k3) = 2f^4P^{ki)P^{k2) + (2 cyclic)]. (94) 



Here, 



M{k) 



2 k^T(k)!,{z) 



(95) 



is the relation in Fourier space between the density and 
the Bardeen potential ^, 



6{li,z) = M{k)(j){li), 



(96) 



where T{k) is the matter transfer function normalized to 
unity as fc — > 0, and g(z) is the linear growth rate of the 
gravitational potential normalized to unity during the 
matter dominated epoch. Further, we define A4i,(fc) = 
M{k)WL{k), Ms{k) = M{k)Ws{k), and so on. Wc can 
then derive the leading contributions in the large-scale 
limit. As shown in App. C, 



X A^z,(A;i)A^i,(fc2)(?!>k(Aki<;^k2) 

= 4/NL<7iC05,L(r), (97) 



where ^4,s,l is the cross-correlation function between the 
density coarse-grained on scale Rl and the Bardeen po- 
tential 6, i.e. 



d?k 
(2^ 



'WL{k) 



d^ki 
(2^ 



(5(k).^(ki)). 



(98) 

In deriving Eq. (97), we have expanded to lowest order in 
k/ki ("squeezed limit" of the bispectrum), with the next 
higher order being suppressed by [k/ki]^ in this limit. 
We will discuss this approximation in Sec. V. 
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Similarly, at leading order in /nl (see App. C), 

X Ms{ki)Ms{k2){(l)\^(l)i,i(t>k2) 
= 2/NLC^5,L(r). (99) 

This result can also be derived by using the well-known 
property of local primordial non-Gaussianity that, in the 
squeezed limit, the local variance of the density field is 
rescaled by <jf>(x), 

(-^^ W) l^(x) = [1 + 4/nl.^(x)] + 0(/^ J, (100) 

and hence 

(y*(x))U(.) =2/nl^(x) + 0(/Sl), (101) 

which immediately leads to Eq. (99). Note that the cor- 
relators involving are independent of the scale -R, for 
local non-Gaussianity, so that the choice of is arbi- 
trary in this case. We will see below how this changes 
for other types of non-Gaussianity. Note that Eq. (99) 
implies that is of order (f), i.e. linear in potential per- 
turbations — in contrast to the naive expectation that 
it is of order 6"^. Finally, we note that (i/»(l)?/*(2)) is 
0(/^l) and hence not included in the following. 



disconnected pieces. In the Gaussian case, Eq. (88) im- 
plies the absence of any connected correlators involving 
Sl and y^ . The powers of then only add zero-lag pieces 
to the previous result Eq. (17), which are absorbed by 
corresponding terms in TV. Thus, the final result Eq. (37) 
does not change in the Gaussian case if we include the 
dependence on y*. 

Defining for convenience 



n,m=0 



we have 



1 



[(/(l)/(2)) - iff 



(105) 



(106) 



where (/) = 0{d^). In the following, we will expand £,h 
to order (5^, and simultaneously to linear order in /nl (as 
long as there are no other sources of non-Gaussianity, 
going to 0(S'^) is also sufficient to retain all terms linear 
in /nl). Through the latter restriction, wo avoid a large 
number of quadratic and higher order terms in y*. We 
have 



and 



1 



C20 



Ur = -^^t + cnC2o{SLy*)<Tl + 0{5'). 



(107) 



(108) 



B. Correlations 

Wc define the estimator for the correlation function 
through Eq. (14) as before. The expectation value then 
becomes 



Hence, £,h becomes 



iUr)) 



1 



" E ;^SSS!<'2(^)^"(^)''2'(2)yr'(2))-i, 

n,m,n' ,m'—0 

(102) 

where '1' and '2' stand for two arbitrary locations sepa- 
rated by a distance r, and we have redefined 



n,m=0 



(103) 



Similarly, we obtain the expectation value of the tracer- 
matter cross-correlation, 

n,m;n-\-m>0 

(104) 

Again, these expressions involve the "bare" bias parame- 
ters Cnm, and the mixed moments of y« which contain 



{ ( (cio6l + coiy, + cnSLy* + ^Sl + ^Sl^^ 
(ciqSl + coiy* -I- CiiSLy* + 



C20 .2 



6 ^/2 



c\q{5l{1)5l{2)) + 2cioCoi(^L(l)y,(2)) 



^20 4 / r \ 2 1 

- -^<^L - CllC2o(<)z,y*)cri j- 

1 

+ cwC2o{5l{1)61{2)) + cwc^ocjI{5l{1)5l{2)) 

+ coiC3oal{5L{l)y*{2)) 

+ 2cnC2o(^L(l)2/*(2))(^L(l)5i(2)) 

+ 2^(5^(1)^^(2))^] +0{6^) 



(109) 



where we have used the symmetry under interchange of 
locations 1 and 2, and 

(2/.(l)5i(2))=3(y.(l)5^(2))ai + 0(/^J 
(5ai)y.(l)5i(2)) = 2(<5^(l)y,(2))(<5^(l)5^(2)) 

+ {5Ly*)al + 0{fl,^). (110) 

Perhaps somewhat surprisingly at first, we have 

to keep these terms whereas terms such as 
(1)2/42)), {5L{l)y.{l)6L{2)) are higher order in 
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/nl and thus dropped. This is simply because the latter 
terms do not have disconnected contributions. 

Note that all completely disconnected terms, i.e. terms 
that asymptote to a constant as r — >■ oo, have canceled as 
expected. We now use the relations derived in Sec. IV A. 
Using Eqs. (97)-(99), we obtain 



1 



+ (2cioCoi + coiC3oo-i + 2cioC2oa-|,) 2fNi,^^s,Lir) 
+ 2cnC2o2hU4>s,L{r)^L{r)] • (111) 

C. Bivariate PBS hias parameters 

In analogy to Sec. II B, we would like to introduce 

a physically motivated bias parameter which quantifies 
the response of the tracer number density to a change in 
the amplitude of small-scale fluctuations, without mak- 
ing reference to any coarse-graining on the scale R^. 
The simplest way to parametrize such a dependence is 
to rcscalc all perturbations by a factor of 1 + e from 
their fiducial value, where e is an infinitesimal parameter. 
For example, for a given realization of initial conditions 
of an N-body simulation, one can obtain a realization 
with a different power spectrum normalization by rescal- 
ing the initial density perturbations by (1 -|- e).^ Clearly, 
the variance of the density field on some scale R, cr^, is 
then rescaled to (1 + e)^(J%- Note that this means that 
the scaled cumulants (5")c/f" are invariant, whereas the 
primordial non-Gaussianity parameter /nl ^ S<j,/P|,, if 
non-zero, scales as (1 + e)^^ under this transformation. 
Specifically, under this rescaling 6l and y* transform as 



y»(x) -)> y,(x) 



(112) 



Note that the parameter in the definition of is just 
a constant normalization, and docs not change vmder the 
e-transformation. This is in analogy to keeping p fixed 
in the iP-transformation in Sec. II B. 

We can then define a set of bivariate PBS bias param- 
eters bNM by generalizing Eq. (31) to 



1 



'JNM = 



D=0,e=0 



(113) 



D=0,e=0 



These parameters can be understood as follows. Given 
infinite volume, the average tracer number density is a 
deterministic function of the mean matter density p and 
the amplitude of the fluctuations (parametrized, e.g.. 



^ Of course, if one initializes using a second-order density field, 
then the second order part needs to be rescaled by (1 -|- e)^. 



through the RM^S of the density fleld on some scale, cr*). 
bNM then denotes the N + M-th joint derivative of this 
function with respect to In J) and (more precisely, s) at 
some flducial values of ~p and a, . Glearly, the parameters 
6jvM are independent of the coarse-graining scale Rl- 

As before, our next task is to derive the relation be- 
tween 6jvM and Cnm- We have from Eq. (92), 



{hu){D,e) = {Fu,Lm £ 



n,m=0 



n!m! 



(114) 



2/* + e 



2 y a? 



We thus have 



(115) 



In particular, 

ho = ^{c,o + fal+0{6l)). (116) 



Further, 



n,m 



^ (coi + C2oai + cMSlV*) + ^ (Si) + 0(6')) . 

(117) 



We can now express the correlation function of tracers at 
this order, Eq. (109), in terms of the PBS bias param- 
eters. In fact, if we are able to reach the analogous re- 
sult to the Gaussian case, i.e. that the tracer correlation 
function is a sum over PBS bias parameters multiplying 
no-zero-lag correlators, we only need to keep terms up 
to order 6^ in bNM, since they always multiply a cor- 
relator of at least order 6"^. Note that when extending 
the treatment to higher order in /nl, it is necessary to 
take into account that j/, transforms nonlinearly with s 
[Eq. (112)]. This means that the bias coefficient multi- 
plying correlators containing say y*(l) will not simply be 
6jV2) but involve a linear combination of 6;vi and 6jV2- 

Let us thus write all mixed "no-zero-lag" terms with 
the appropriate bNM in front, at order 5^, /nl- We obtain 



a(r-) 



^ bUUr) + ^Ci(r) + 2&io6oi(^l(1)2/*(2)) 



&2o6ii(<5L(l)y»(l)'5i(2))r 



Here 
der. 



+c?(<5^/^L)• 

we have used the fact 
{5l{1)SI{2)), {6l{l)y42)), 



that at this 
{Sl{1)51{2)), 



(118) 

or- 
and 



(y*(l)(5|^(2)) have no no-zero-lag pieces. Note also that 
('5L(l)y*(2)) = (i5L(l)2/*(2))nzi- Plugging in the expres- 
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sions for 6jvm at the relevant order, we obtain 

Ur) = ^{ (c?o + CioC30^i) ^L{r) + ^^1(0 

+ 2 (^ciocoi + ^coicrl + cioC2oal^ 2/nl^<#.5,l(?') 

(119) 

We easily see that this agrees identically with Eq. (111). 

Thus, the bivariatc PBS parameters which wc have de- 
fined in a coarse graining scale-independent way absorb 
all coarse graining-scale dependent terms in the "bare" 
bias parameter expansion Eq. (102), in particular the 
term ciC2((5l(1)<5|(2)). We expect this to hold to any 
order in the bare bias parameter expansion, although a 
proof is beyond the scope of this paper. 

Thus, the introduction of the bivariate bias parameters 
Eq. (113) and the resulting expression Eq. (118) achieved 
exactly what we had wanted. In particular, the leading 
effect of local primordial non-Gaussianity is quantified 
by 6oi, the response of the mean number density of trac- 
ers to a rescaling of the amplitude of initial fiuctuations. 
The term ciC2{Sl{1)S^{2)) on the other hand is seen as 
an artifact of the bare bias expansion which is absorbed 
in the renormalized parameter boi. Apart from the clear 
physical interpretation, this reordering of the perturba- 
tive expansion is also manifestly convergent: higher or- 
der terms are guaranteed to be suppressed by powers of 
^L{f) and /nlCi;)<5.l(''), which only need to be small on 
the scale of observation r for the perturbative expansion 
to be valid. 

This also remedies a worrying issue with the local 
bias expansion in the non-Gaussian case: evaluation of 
Eq. (36) shows that higher order terms ("loop correc- 
tions") become comparable to or larger than the lead- 
ing order expression &i^L(r) on sufficiently large scales, 
which would indicate a breakdown of the perturbative 
expansion. The bivariate expansion on the other hand 
leads to an expansion in which higher order terms are 
consistently suppressed [Eq. (118)], i.e. all dominating 
terms are actually lowest order ("tree- level"). For suf- 
ficiently large values of /nl, one might need to include 
higher order terms in that parameter. Nevertheless, the 
expansion will remain convergent. 



D. Non-local non-Gaussianity 

Wc now consider the generalization of the results of 
the last section to arbitrary quadratic non-Gaussianity, 
i.e. non-Gaussianity that is described to leading order 
by a 3-point function. The correlators that are relevant 
for the tracer two-point correlation in the non-Gaussian 
case, Eqs. (97)-(99), are determined by the behavior of 
the bispectrum in the squeezed limit, corresponding to 
triangle configurations where one side is much smaller 



than the other two. For scale-invariant bispectra, we can 
write the bispectrum in this limit as 

(120) 

with A, a being constants (more general shapes can be 
constructed by linear superposition of bispectra with dif- 
ferent Aj, (li). Local, folded, and equilateral shapes cor- 
respond to a = 0, 1, and 2, respectively. Eqs. (97)-(99) 
then generalize to 



(5. 



{5L{i)y*{2)) 



X Ml{kr)Ak'^P4,{k)k:['^P^{k{) 
= Aa\^LU.s,L{r) (121) 

= 2\''^W^/ 



X Ml{ki)Ak"'P4,{k)k^"P^{kx 



(122) 



where we have defined the general spectral moment 



'n,X 



:k^P{k)\Wxm, (123) 



and the correlation between a non-local function of (p and 
the density field. 



(Pk 
(27r)3 



k''M{k)P^{k)WL{k). (124) 



Again, Eqs. (121)-(122) are valid at leading order in 
the squeezed limit {k <C k\, with corrections going as 
{k/kiY). Inserting these expressions into Eq. (109), and 
using Eq. (110) we obtain up to 0{5^, /^l) 



a(r-) = 



^2 



(c?o + cioC3oai)a(r) + |^a(r)^ 



(^2cioCoi^ 



■ C01C30- 



Oi,s 2 I 2 



2a? 



+ 2ciiC2o^^^^^„5,i(r)^i(r) 



(125) 



Inspection shows that the bivariate PBS parameters de- 
fined in Sec. IV C cannot absorb the i?L-dependent term 
from (5^(1)^1(2)). This goes back to the fact that in 
the presence of a bispectrum of the form Eq. (120), the 
small-scale perturbations are not rescaled uniformly, but 
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rather in a scalc-dcpcndcnt way: the squeezed-hmit re- 
sult Eq. (100) gcncrahzcs to 



where 



As in the case of local non-Gaussianity, we now write 
all mixed "no-zero-lag" correlators with the appropriate 
h^M in front, up to 0{5^ , f"^-^). Due to the factor of 



('^s(^))|0^(x) ~ ^« + ^f-a.s^aW' (126) ^^^2 -j^ ^j^g transformation of under the scale- 



.« = / 



(2^ 



A;"«!)(k)e 



ikx 



and hence 



(127) 



dependent rescaling Eq. (129) (at lowest order), we have 
to divide by that factor when multiplying correlators in- 
volving y*. We obtain 

ih{r) = bUdr) + ^^i(r-) + 26io&oi^ (5^(1)2/* (2)) 

— r^ a 



(128) 



Thus, the transformation of the density field following 

Eq. (112) is not the relevant one any more. Instead, we 
need to rescale the density field through 

(5(k) (l-h£fc-")5(k), (129) 

so that 

^l(x) 6l{x) + £ ^_a,L(x) 



where, in analogy to Eq. (127) 



(130) 



(131) 



Note that {5x5-a,x) = <y-c,x^ and {S_^ x) = '^-2a,x- 
We will continue to assume that the tracer density de- 
pends on the small-scale density field only through the 
variance on some scale i?*, parametrized through y*. We 
again define b^M through Eq. (113), but with the trans- 
formation Eq. (129), so that these bivariate bias param- 
eters will in general be different from those in Sec. IV C. 
As before, our next task is to derive the relation between 
bNM and c„m- We have from Eq. (92), 



n\m\ 



{dL + s6-oc,L)SL+D 



s s 
y* + -^Ss6_a,s + ^<5- 



We obtain 



1 C„,; 

\f ^-^ n'rr 



M '■ — ' n\m\ 



m 



coi- 



+ C20T±a,L + O{6',f^^) . (132) 



-I- b2ob 



-(5^(l)y,(l)5i(2))„,, 



= bUL{r) 



^^eLir)+bioboiA^<p„sAr) 



b2obiiA^^^s,L{r)^L{r). 



(133) 



Note that the final result is explicitly independent of 
the scale Rl (as long as r is sufficiently large so that 
the smoothing effect on £,4,^s,L{r) is negligible), 

whereas ((5i,(l)y*(2)) itself is not since it depends on the 
spectral moment cr^ [Eq. (122)], which in turn depends 
on al [Eq. (89)]. 

Inserting the expressions for 6jvm at the relevant order, 
and using Eqs. (121)-(122), we have 

a(r) = ^1 {clo + ciocsoai) Ur) + f^^i(r) 



2cioCoi + 2cioC2oCr^a r- 



X ((5i(l)j/.(2)) 

+ C20Cll(^L(l)y*(l)<5i(2))nzl| 

= ^{(^?o + CioC3oai)eL(r) + |'ei(r) 

+ [2cioCoi -I- C30C01(t|] -^^A.^4,„s,R{r) 

+ CioC2oO'-Q,L^^<i>„<S,-R(r) 

+ 2c2oCii^L{r)^^A^4,^sAr)y (134) 

This agrees exactly with Eq. (125). The key differ- 
ence of the expansion of ^/^ in terms of renormalized 
bias parameters in the case of non-local primordial non- 
Gaussianity, Eq. (133), from the corresponding result for 
local non-Gaussianity Eq. (118) is that the bivariate bias 
parameters are now defined with respect to the scale- 
dependent rescaling of the density field, Eq. (129), rather 
than a scale-independent rescaling. We find that it is 
sufficient in the large-scale limit, even in the case of a 
non-Gaussianity of general shape, to describe the coarse- 
grained tracer abundance as a function of (5l(x) and 
?/*(x) in order to absorb the dependence on the coarse- 
graining scale i?i into the bivariate PBS bias parame- 
ters. However, the actual definition of the renormalized 
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bias parameters depends on the shape of primordial non- 
Gaussianity, in particular the scaling with ki/ks in the 
squeezed limit. 

We can thus summarize our findings regarding the ef- 
fect of a primordial bispectrum on the two-point corre- 
lations of tracers (non-Gaussian scale-dependent bias) as 
follows: 

• For local primordial non-Gaussianity, it is sufficient 
to include the dependence of the tracer density hh 
on the local amplitude of small-scale fluctuations 
5s through the variance on some scale ij, . Further- 
more, the scale i?* (and whether the dependence on 
5s is actually through the variance on several scales) 
is irrelevant, as all perturbations 5s are rescaled 
uniformly. 

• For non-local separable bispectra as in Eq. (120), 
it is still sufficient to parametrize the dependence 
of hh on the amplitude of small-scale fluctuations 
through the variance on a single scale . However, 
the value of the scale i?« now matters as is modu- 
lated by an amount that depends on [Eq. (128)]. 
In particular, if the tracer number density were to 
depend on the variance of 5g on several different 
scales, then the PBS bias parameter 6oi will be a 
linear combination of these different dependencies 
with relative weights controlled by a, i.e. the shape 
of the bispectrum. 

• For non-separable bispectra, the renormalization 
approach we describe here is not able to remove 

the i?i-dcpendence in the tracer correlation fimc- 
tion. However, such shapes can typically be well 
approximated by a linear superposition of separa- 
ble shapes (see e.g. [43]), which then allows the 
renormalization to proceed as described here. 

Thus, we find that in general, a given tracer will re- 
spond differently to different shapes of primordial non- 
Gaussianity, i.e. feoi (and &7vm hi general with M > 0) 
depends on the tracer as well as the shape of the primor- 
dial bispectrum. In the following we will study this in 
the context of simplified models of tracers. 



E. Universal mciss functions 

We begin with a generalization of the universal mass 
function discussed in Sec. II C. We write the mean abun- 
dance of tracers as 



to lowest order as 



nh =nh [p, (^*,J*) 



(135) 



where the Jacobian J* is defined in Eq. (39). That is, 
fih is given as a function of the mean density of the Uni- 
verse and the variance of the density field smoothed on 
a scale i?*, as well as its derivative with respect to scale. 
Under the generalized rescaling Eq. (129), a* transforms 



l + £- 



(136) 



while the Jacobian transforms as (see also [34]) 



l + 2e- 



'dlncr^^ _ 
rflnii. 



din a? 



dlncri^^ 
rflncrj 




(137) 



Here we have used d/d\nR^ = 2 J* d/dlna^. Note that 
in the local case where a = 0, the local Jacobian is not 
affected by long-wavelength modes. Using Eq. (135), we 
can then derive the leading non-Gaussian bias through 
Eq. (113): 



boi = 



duh d In cr* 




duh d In J* 
d In J* de 



I din a? 



dlncrj 



(138) 



Here, 6oi(q^ — 0) is the PBS bias parameter quantify- 
ing the effect of local primordial non-Gaussianity for a 
tracer following Eq. (135), and we have assumed that 
the tracer density scales linearly with the Jacobian as 
expected physically. For such tracers, the bias parame- 
ters quantifying the response to general non-local non- 
Gaussianity (in the squeezed limit) are thus directly re- 
lated to those for local non-Gaussianity. In particular, 
we recover the results of [34], who first pointed out the 
contribution by the Jacobian J* . 

We now specialize Eq. (135) to a "truly" universal mass 
function [Eq. (38)], 



nh= P f{i^c) J* , 



(139) 



where f{vc) is in general an arbitrary function of Vc- The 
results relating 601(0^) to feoi(a = 0) of course also hold 
in this case. However, the specific form Eq. (139) further 
allows us to connect 601 (o^ = 0) to the linear PBS density 
bias: 



310 



6oi(a = 0) = 



1 dfih 1 df 

hh 9 In p (T* dvc 

1 dhh _ 5c df 
hh d In (7* cr* duc 



= 5J)w. (140) 



Note that here 5io is the Lagrangian bias, which is why 
we have not included the derivative with respect to In p of 
the p prefactor in Eq. (139) (see also Sec. II C); again, the 
effect on J* vanishes for a = 0. This is the original rela- 
tion between the density bias parameter and the response 
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to primordial non-Gaussianity derived in [13, 14, 41]. Wc 
point out that these resuhs differ from those of [44] , who 
considered the effect of primordial non-Gaussianity on 
tracers with local Lagrangian biasing. There, the en- 
tire leading order effect of primordial non-Gaussianity 
is encoded in a scale-dependent C2(ki,k2). Thus, a 
parametrization of C2(ki, down to very small scales is 
necessary in order to predict the amplitude of the scale- 
dependent bias. This is in contrast to the approach pre- 
sented here, where one introduces a local dependence on 
the small-scale fluctuations which absorbs the term pro- 
portional to C2 into a renormalizcd &oi, which is a single 
number. As a result, the prediction of [44] yields a depar- 
ture from Eq. (140) for universal mass functions which 
depends on the precise form of /(^c)- While Eq. (140) 
has been both supported [13, 41, 45] and disputed [46] by 
simulation resiilts, these different predictions are clearly 
resolvable with sufSciently large simulations. In partic- 
ular, our prediction for the scale-dependent bias for a 
general tracer. 



which is independent of any assumptions on the mass 
function of the tracer, provides a rigorous test of our 
approach which can be applied to simulations. 

V. SUMMARY &£ DISCUSSION 

We have shown that the expression of tracer correla- 
tions in terms of i?L-indepcndcnt renormalizcd bias pa- 
rameters 6jv absorbs all zero-lag correlators present in 
the expansion of the tracer correlation function in terms 
of the bare "scattcr-plot" bias parameters c„. We have 
shown this to all orders for an arbitrary density field. 
While the proof only applies directly for the auto- and 
cross-correlation functions and pure density biasing, we 
expect the result to hold in the case of higher A''-point 
functions and multivariate biasing as well (analogously to 
the resummed multipoint propagators of [10, 11]). Our 
key result is a rigorous definition of the renormalizcd bias 
parameters in terms of derivatives of the mean number 
density of tracers with respect to the background density 
(we call these "peak-background split" bias parameters 
since their definition is closely related to the commonly 
adopted definition of PBS biases [30, 47, 48]). It is im- 
portant to stress that this exact definition is entirely in- 
dependent of the nature of the tracer considered. There- 
fore, it provides a rigorous framework in which further 
assumptions for or modeling of the bias parameters, for 
example from the excursion set, peak model, or halo oc- 
cupation distribution, can be embedded. 

Our rcsiilts go beyond previous work on renormalizcd 
bias parameters [8] in two ways: first, we show that our 
result is valid to all orders; second, we rigorously con- 
nect the renormalizcd bias parameters with the peak- 
background split. We further show that the renormal- 



izcd bias parameters in the tracer auto-corrclation and 
the tracer-matter cross-correlation agree to all orders. 
We also expect this to be the case for higher iV-point 
functions, although this remains to be shown. We can 
summarize this reasoning as in line (a) of Tab. I: the ex- 
pression of tracer correlations in terms of the bare biases 
c„ is i?i-dependent at each order due to disconnected 
correlators (for example c|cr|). This i?_L-dependence is 
then resummed into _Ri-independcnt bias parameters bjsf 
which are defined with respect to a uniform increase in 
the matter density. 

The underlying assumption in this result is that the 
clustering of tracers is entirely determined by their de- 
pendence on the local matter density. This is not ex- 
pected to be a good assumption in general. However, 
our result provides another invaluable tool: whenever the 
renormalized expression in terms of no-zero-lag correla- 
tors exhibits a residual dependence on Rl, we conclude 
that a biasing purely in terms of matter density is not 
sufficient. 

We first encounter this in the case of the smoothed 
matter correlation function ^L(r), which depends on Rl 
if ^(r) has structure on scales smaller than Rl (line (b) 
in Tab. I). In this case, we are led to introduce bias pa- 
rameters with respect to the curvature (Laplacian) of the 
matter density field. In Fourier space, this corresponds 
to a scale-dependent bias oc k^. If we further include 
bias parameters with respect to higher derivatives of the 
density field, we can in fact entirely absorb the effect of 
smoothing on ^(r) [App. B]. Again, this is regardless 
of the nature of the tracer and the shape of the matter 
correlation function. Of course, for a smooth correlation 
function, it is usually sufficient to keep only terms involv- 
ing the lowest few derivatives of the density field. The 
renormalized biases with respect to the curvature corre- 
spond to derivatives of the mean tracer abundance with 
respect to a constant shift in the curvature of the density 
field (Tab. I). As an example, these bias parameters are 
easily derived for peaks of a Gaussian density field from 
the results of [39]. We show that the bias parameters 
obtained in this way indeed match the scale-dependent 
biases derived in the full, direct calculation of peak corre- 
lations [36]. In this context, it is important to point out 
that the curing of i?L-dependencies, such as that from a 
smoothing of the correlation function, is a sufficient con- 
dition for having to introduce an additional dependence 





f?i,-dependence 


local 
quantity 


Transformation defining 
PBS bias parameter 


(a) 




5l 


p + D'p 


(b) 








(c) 


ciC2{c5L(l)<5i(2)) 


y* 


5^ {l + e)5 



TABLE I: Summary of renormalization procedures intro- 
duced to remove various dependencies of tracer correlations 
on the coarse-graining scale Rl- The variable y* is defined in 
Eq. (90). 
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of the tracer density on properties of the matter density 
field. However, it is not a necessary condition — specific 
tracers might also exhibit additional dependencies not re- 
quired by renormalization. One example is peaks of the 
matter density field, which also exhibit a dependence on 
quantities such as (V(5)^ [40]. Of course, it is straight- 
forward to include these additional dependencies in the 
formalism described here, by defining renormalized PBS 
bias parameters through suitable transformations of the 
density field. 

In the case of a non-Gaussian density field, we find that 
the tracer correlation function for pure density biasing 
acquires a strong dependence on Rl if long- wavelength 
modes are coupled to short wavelength modes (line (c) 
in Tab. I). The most well known example of this kind 
is primordial non-Gaussianity of the local type [18]. In 
this case, we have to add a bias parameter with re- 
spect to the amplitude (variance) of small-scale fluc- 
tuations. The renormalization procedure then absorbs 
the i?L-dependent terms such as CiC2(5i(l)(5|(2)), and 
the resulting bivariate bias parameters are given by the 
derivatives of the mean tracer density with respect to 
the background density and (essentially) the amplitude of 
the initial power spectrum — both clearly i?z,-independent 
quantities. Effectively, we obtain an expansion closely re- 
lated to that of [41], although we did not need to drop 
any terms or make approximations beyond the large-scale 
limit (which allows us to evaluate the bispectrum in the 
squeezed limit). 

We also generalize the results to any form of primordial 
non-Gaussianity given through a bispectrum of potential 
perturbations. In fact, this provides a good example for 
how this renormalization approach pays oS: we obtain a 
fully general and exact result (in the large-scale limit), 
in which the renormalized scale-dependent bias parame- 
ter depends on the precise shape of the non-Gaussianity 
as well as the nature of the tracer. Assuming that the 
tracer abundance only depends on the variance of the 
small-scale density field on a single scale, we can how- 
ever relate the scale-dependent bias parameter for an ar- 
bitrary general shape to that for local non-Gaussianity. 
Further, we can be more restrictive and assume that the 
tracer follows a universal mass function. In this case, we 
can relate the scale-dependent bias parameter to the bias 
parameter with respect to density (as in [13, 14, 41]). 

The general procedure also carries over to primordial 
non-Gaussianity described by higher iV-point functions. 
For example, a non-zero trispectrum which couples long- 
to short-wavelength modes will introduce a significant 
i?L-dcpendence in the tracer correlation function through 
the term cic^{Sl{1)S^{2)) . In order to remedy this, we 
need to explicitly take into account the dependence of the 
tracer density on the local skewness (S^) of the density 
field, which then yields a corresponding scale-dependent 
bias contribution (as shown in [34, 49]) which absorbs the 
iii-dependent terms. 

The main caveat to our results is that we have worked 
in Lagrangian space throughout. While we expect that 



the general approach will also be applicable in Eulerian 
space, the effect of gravitational evolution will in general 
introduce several further dependencies of the tracer den- 
sity on the environment, for example velocity and tidal 
field biases [19-21]. We leave this for future work. 

We have also only considered observables in real space. 

The main obstacle in transforming to Fourier space is 
the issue of stochasticity in the tracer density field and 
its scale dependence, which contributes to correlations 
at all k in Fourier space (although the contributions will 
asymptote to a constant in the low- A; limit). Thus, a well- 
defined model for correlations on small scales is a neces- 
sary prerequisite for a rigorous understanding of Fourier- 
space correlations. 

Further, we have restricted the treatment here to two- 
point correlations of tracers. The main reason for this 
is simplicity; we expect no major obstacles in generaliz- 
ing the results to higher TV-point functions in Lagrangian 
space, such as the tracer bispectrum with non-Gaussian 
initial conditions. In order for this to be useful however, 
non-Gaussianities from gravitational evolution will also 
have to be included [50-55]. 

Finally, in the case of primordial non-Gaussianity, we 

have only considered linear terms in /nl, and restricted 
to the large-scale limit where the bispectrum is evaluated 
at lowest order in the "squeezed-limit" expansion. The 
extension to higher powers of /nl is straightforward. The 
second approximation captures the main effects on large 
scales, since the subleading term is suppressed by kf/kg, 
where ki ~ 1/r is the scale on which we measure correla- 
tions, and ks ~ 1/i?* corresponds to the small-scale fluc- 
tuations. For example, in case of local non-Gaussianity, 
the subleading term is expected to lead to a small ap- 
proximately scale-independent bias. 

These caveats notwithstanding, we hope these results 
provide the starting point for a rigorous treatment of 
biasing of general tracers in the context of cosmological 
perturbation theory. 
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Appendix A: Proof of Eq. (36) for a general 
density field 

Let us denote as n„ the set of all partitions of the set 
{ 1, 1, 1 }, (Al) 

n times 

where the elements of the set are considered distinguish- 
able. We call the elements B of any given tt € n„ 
"blocks", with tt having |7r| blocks where |7r| is the cardi- 
nality, or number of elements, of tt. Clearly, |7r| < n, and 
the blocks of any partition in n„ satisfy 

For example for n = 4 there are four distinct partitions 
with one block Bi with = 1 and one block B2 with 
1-62! = 3. Then, the moment for an arbitrary density 
field 6 is given in terms of the cumulants (connected cor- 
relators) by 

r>= E n<'^"")- (A3) 

7ren„ BeTT 

For example, in this sum the trivial partition tt = 
{{1, 1, 1}} (with a single block B with \B\ = n) 
corresponds to (i5")c- Note that since {S) = 0, any par- 
tition where \B\ = 1 for any B G tt yields a vanishing 
contribution. 

Similarly, let us denote asUnm the set of all partitions 

of 

{1, 1, 1, 2, 2, ...,2}. (A4) 

^ V ' V ' 

n times m times 

We can then write 

(W)= E IKII'^"), 

= E n(c^''^<^2^^''')c, (A5) 

pen„,„ Bep 

where in the first line a runs over the elements of the 
block B, and in the second line we have defined as ni{B) 
the number of elements '1' in block B, and correspond- 
ingly for n2{B) (so that ni{B) + n2{B) = \B\). This 
simplifies the result since the cumulants are independent 
of the order of products of 5i and ^2, and only depend 
on the overall power of each. 

Our goal is to reorder the sum in Eq. (A5). We assign 
two further numbers (non-negative integers to be precise) 



to each p e 

„2(B)=0 rii(B)=0 

N,ip)= n,{By, N2{p)= ^ n2{B). 

Bep Bep 

(A6) 

In other words, for a given partition p, Ni (p) counts the 
number of elements '1' that are in blocks that only con- 
tain '1', while N2{p) counts the number of '2's that are in 
blocks only containing '2'. These numbers are of course 
uniquely defined for each p. Moreover, TVi, N2 define a 
partition of Iln.m, i.e. each p G Iln.m is member of one 
and only one subset of Tin,m defined as containing all p 
with a specific value of Ni and N2. Equivalently, the 
relation p a defined for any p,a G Hn,m through 

p^a ^ Ni{p) = Ni{a) A Niip) = Mcr) (A7) 

is an equivalence relation on n„_„. We can then split the 
sum in Eq. (A5) into sums over these disjoint subsets of 

n m Ni(p) = k- N2{p)=l 

fe=o ;=o peUn m Bep 

(A8) 

Consider the sum over all partitions in one of these sub- 
sets: 

Ni(p)=k; N2{p)=l 

pen„,„ Bep 

A partition p E Iln.m can be thought of as one specific 
way of distributing n black balls and m red balls into ar- 
bitrarily many (initially empty) boxes. These boxes cor- 
respond to the cumulants in Eq. (A5) (of course empty 
boxes are trivial, because they yield 1 in the product in 
Eq. (A5); boxes with only one ball lead to a zero con- 
tribution). The sum in Eq. (A9) runs over all possible 
ways of distributing these balls that have exactly k black 
balls which are in boxes with only black balls, and I red 
balls which are in boxes with only red balls. Correspond- 
ingly, the remaining n — k black and m — I red balls are 
in boxes with both black and red balls. There are (n k) 
ways of selecting k black balls out of n, and (m I) ways 
for the red balls. Given such a selection of k out of n 
and I out of m, the sum in Eq. (A9) thus runs over all 
ways of partitioning k black balls into boxes, m red balls 
into a different set of boxes, and finally n — k black and 
m — / red balls into a third set of boxes such that each 
of these boxes contains at least one black and one red. 
Mathematically, we can write Eq. (A9) as 
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JVi(rt=fc;iV2(p)= 



Bep 



n 



E n i^rh 

j2 n^c^^^c^^^) 



m 



^2 611; B2e7r2 



"l(S)rn2(B) 
"2 



(AlO) 



Here, we have used Eq. (A3), and defined the subset Ilf^l^ 
of the set of all partitions Iin,m through 

Ul% = {pe n„,„ : 7Vi(p) = = N2{p)} 

= {pen„,„: VB e p ni(B) > An2(B) > 0|. 

In other words, H"^^ contains all partitions in which each 

block has at least one element '1' and at least one element 
'2'. That is, for any p £ njj^)^, the product of correlators 



Bep 



(All) 



does not contain any zero-lag pieces. It is then natural 
to define the no-zero-lag correlator for a general density 
field as 



penn^L, Bep 



«i(B)r"2(B) 



)c , (A12) 



which reduces to Eq. (25) for a Gaussian density field. 

Finally, we can sum Eq. (AlO) over all values of 
^i(p), N2{p) to obtain 

{SiST) = EE ) {Si) (7) {Sl){6r'ST-\.i ■ 



fe=0 1=0 



(A13) 

We are now ready to prove Eq. (36) for general non- 
Gaussian matter density fields. Plugging the relation be- 
tween &JV and c„, Eq. (33) into Eq. (36), and relabeling 
N^n-N,M^m-M, yields 

n— 1 m— 1 



iUr))=^ E SSE E 



nlmi 




X )(52-^(l)<5r-^^^(2)>„.i 



^ oo 

= 7V^ E 

n,m— 

oo 

= Jf2 E 



n,m=0 



^(<5r(l)<5r(2))-l. (A14) 



In the second line, we have used Eq. (A13) and subtracted 
out the completely disconnected contribution 
(which is not included in Eq. (36) as we start the sum 
from = 1, M = 1). In the third line, we have used 
Eq. (22). This is identical to Eq. (17), and thus proves 
the relation Eq. (36) between PBS bias parameters and 
two-point correlations for a general non-Gaussian density 
field at all orders. 



Appendix B: Curvature hias to higher order 

Let us consider a general spherically symmetric filter 
function, 



(Bl) 



where /(y) is a dimensionless function defined on [0,oo). 
We have pulled out a factor of (47ri?|^)~^ for convenience. 
The Fourier transform of Wl (x) can then be written as 



n=0 



Jo 



(B2) 



The normalization constraint / rf^xWj:,(x) = 1 is equiv- 
alent to /o = 1. Since i?^ is an arbitrary parameter, we 
can further choose one of the /„ with n > 1 to assume 
some desired value. Specifically, in Sec. Ill we have cho- 
sen Rl so that /i = 6 and hence W{k) = 1 — k'^R\ + • • • . 

The smoothed correlation function can then be written 
exactly as 

oo f f p2(n-t-m) 



n.m=0 



Note that the factor of (—1)"+™ from the expansion of 
Wl cancels with the ^^(n-i-m) fjom converting powers of 
k into derivatives. 

We can generalize the transformation Eq. (61) as fol- 
lows: 

oo 

<5a(y) = %) + aE|^|yr, (B4) 



JV=0 
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OO ^ 

^L,a(0) =^L(0)+a^ ( 



(B5) 



where a and g^r are dimensionless parameters, and the 
case A'^ = corresponds to the transformation used to 
derive the PBS density bias (with — \ and a — D, 
Sec. II B). Using the expansion of the filter function, we 
obtain 

2 A' 

J 5iv/jv- 

Further, using that 

V2|y|2^ = (27V+l)2iV|y|2^-2 

(2jV + l)! 
(27V- 2n+ 1)!'-" 
for n < N, wc obtain in analogy to Eq. (B3) 
V2"<Si,„(0) = V2",5i(0) 

" 'i?iy(^""'_(27V+l)! 



y2n| |2iV 



|2(JV-n) 



(B6) 



a \ ^ 



N= 



{2N-2n + iy. 



QnIn-u- (B7) 



The transformation Eq. (63) corresponds to gi 
with all other g^ equal to zero. Hence, 

V25i,„(0)=V2<5i(0) + |3!^ 



= V25l(0) 



a 



1/6 



(B8) 



as intended (with all higher derivatives being unaffected). 
More generally, we can define gN = 1/(2A''+1)! for a fixed 
N , with all other g„ = 0, so that 



V'"<5i.„(«,(0) = V^"5i(0) 



+ 



1 



\ I J (27V-2n+l)! 
where n < N. In particular. 



N—ni 



(B9) 



^L,aw(0)=<5L(0) + aW 



-.1 



N 



V^^^L,„w(0)=V^^540) + 
Now, if we write 



I 

(2N ■ 



n,(x)=F,x({v2"5L(x)}„^,;x 

we can define the generalized, bare linear curvature bias 
parameters (at linear order in V^"(5) through 



(2iV+ 1)!- 

(BIO) 

(Ell) 



_ 1 dFu,L ({V^"^L(x)}::^3;x) 

CV2"5 = 



Further, we can define PBS bias parameters as 



a(")=0 



(B12) 



(B13) 



N 



n=0 



N 



^Cv2 
n=0 



.2(N-n) 



R 



'*(2Af-2n + l)!-^'^-"' 



(B14) 



where we have taken out powers of £ to make the expres- 
sion for the correlation function below simpler. Note that 
the usual PBS density bias is obtained as special case for 
N = 0. 

In terms of the bare bias parameters, and to linear 
order in matter correlations, we can write the tracer cor- 
relation function as 

OO 

a(r-)= J2 Cv-5Cv2^5V2("+™)a(r) 

n.m— 



E 

n.m— 



E 



f f p2(JV+M) 

JnJm^l 

(27V + 1)!(2M + 1)! 

JV,M=0 ^ ^ ^ ^ 

OO N,M 

N,M=On,m=0 

r r p2(N-n+M-m) 

^ {2N -2n + 1)!(2M - 2m + 1)! 



y2(JV+M+n+m)^^^^ 



(B15) 



y2(iV+M)^(^-)_ 



In the second line, we have shifted the sum by defining 
N = N + n, M = M + m, while in the third line we have 
reordered the sum over n, m and A'^, M. 

Now let us write down the expected result in terms of 

PBS bias parameters, i.e. assuming that the parameters 
by2Ng have absorbed all dependencies on Rl'- 

OO 

a(r-)= E &v-5&v2MiV2(^+^)^(r), (B16) 

Af,M=0 

where we have again restricted to linear order in matter 
correlations. Using Eq. (B14), this yields 



OO N,M 

E E Cv2"5Cv2m5 

N,M=0 n,m=0 



R 



2{N+M-n-m) 



fN-nfm- 



{2N -2n+ 1)!(2M - 2m + 1)! 



V2(^+^)^(r), 



which agrees with the exact result Eq. (B15). Thus, the 
PBS bias parameters defined in Eq. (B13) are indeed 
able to completely absorb the effects of smoothing on 
the correlation function (at linear order). Moreover, if 
we assume that there is a characteristic scale £ describing 
the dependence of the tracer density on the derivatives 
of the density field, i.e. 

b^2Ms ~ (B17) 

then the expression Eq. (B16) is an expansion in terms 
of 



(B18) 



If the matter correlation function does not have signif- 
icant structure on scales below £, then this quantity is 
progressively suppressed at higher N and Eq. (B16) is 
rapidly convergent. 
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Appendix C: Derivation of Eqs. (97)-(99) 



while 



In this appendix we derive the squeezed-hmit expres- 
sions Eqs. (97)-(99). We begin with the correlator 
Eq. (97): 

(«lMi(.)>^/|^e-M.,.,/i|/|| 

Using the definition of the bispectruni, 

(<^k<^ki0k.) = (27r)3jB(k + ki+k2)B^(k,ki,k2), (CI) 
we have 



{5,m{2)) = J j^.e^'^-Mdk) J ^,Mdk. 

X 7WL(|k + ki|)5^(fc,A;i,|k + ki|). 

(C2) 

We now expand the integrand in powers oi q = k/ki. 
Further, we define /x = k • ki. We then obtain for any 
function F{k) 

i^(|k + kil) = F(fci) [1 + {2qn + q^)a + 2gV'&] + 0(9'), 
^ d\o.F{ki) 
d\nk\ 



d'^lnF{ki) ainF(fci) 



B4,{k, ku |k + ki) = 2/NL{P0(fc) [P^{ki) + P^dk + ki|; 

+ P0(A:i)P4|k + ki|)} 
= 2fj,Mk)P^{ki) [2 + {2qn + [1 - 2/i2]g2)„^] + ^(^3)^ 

where = Ug — 4 and we have assumed a pure power- 
law P,p{k) for simplicity. Inserting these expressions into 
Eq. (C2), we see that the terms oc q/j, vanish once the 
integral over fi is performed. Thus, we obtain 



(5i(l)5i(2))=4/NL / ^e''^-''ML{k)P4k) 



d{lnkif dlnki ' 



(C3) 



d^k 
(27r)3 

^MlikMk^) [1 + Oiq')]. 

The ki integral yields a£, leading to Eq. (97) with cor- 
rections suppressed in the large-scale limit by {k/ki}^ ~ 
{kRLf. We now turn to Eq. (99). Since 

(54l)y.(2)) = ^(5^5,2(2)), (C4) 

this just differs from Eq. (97) through the prefactor and 
the different filter function, and the above results imme- 
diately lead to the third line of Eq. (99). 
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